Dispersive magnetic excitations in the S = 1 antiferromagnet Ba3Mn2 08 
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We present powder inelastic neutron scattering measurements of tlie S — 1 dimerized antiferro- 
magnet Ba3Mn208. The T — 1.4 K magnetic spectrum exhibits a spin-gap of A « 1.0 meV and 
a dispersive spectrum with a bandwidth of approximately 1.5 meV. Comparison to coupled dimer 
models describe the dispersion and scattering intensity accurately and determine the exchange con- 
stants in Ba3Mn208. The wave vector dependent scattering intensity confirms the proposed 5=1 
dimer bond. Temperature dependent measurements of the magnetic excitations indicate the pres- 
ence of both singlet-triplet and thermally activated triplet-quintet excitations. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Et 



I. INTRODUCTION 

Low-dimensional and gapped quantum magnets based 
upon strongly coupled spin pairs or dimcrs with weaker 
interdimer interactions have become especially topical 
systems. This is primarily due to the relevance of exper- 
imentally accessible quantum critical pointsii^i^ii^. For 
antiferromagnetic intradimer exchange, the ground state 
of such systems is a product of singlets, but strong mag- 
netic fields can close the spin-gap to excited triplet states 
via Zeeman splitting of the triplel^i^ii. Such systems 
thus provide an elegant realization of a lattice gas of 
hardcore bosons in which the external magnetic field 
plays the role of the chemical potential and the inter- 
dimer coupling determines both the kinetic and poten- 
tial energy of the delocalized triplets^i^. Depending on 
the balance of these energy scales the triplets will ei- 
ther crystallize or condense at low temperatures^^ or, 
under the right set of conditions, form a supersoliclii. 
Ba3Mn208 is a particularly promising candidate mate- 
rial for the detailed study of magnetic field dependent 
quantum critical points. Ba3Mn208 has been identi- 
fied as a 5 = 1 dimerized antiferromagnet and the low- 
temperature phase diagram has been examined using 
thermodynamic measurement a^^'^'^i*^ . However, a mea- 
sure of the dominant exchange constants has been no- 
tably absent. 

Low-field magnetic susceptibility and pulsed field 
magnetization measurements of Ba3Mn2 08 agree with 
weakly coupled S = 1 dimer models with antiferromag- 
netic dimer exchange Jo between 1.50 and 1.65 meV 
and a zero-field spin-gap of A = 0.97 meV— i^^. Un- 
fortunately, comparisons to isolated dimer models only 
yield results in terms of weighted sums of interdimer ex- 
change. Specific heat measurements on powder samples 
have been performed in magnetic fields up to fXoH = 29 T 
to describe the low-temperature magnetic phase dia- 
gram of BaaMnzOgi^. At T w 0.6 K, a critical field 
of fJ-oH w 11 T is required to induce magnetic order- 




FIG. 1: (Color online) Crystal structure of Ba3Mn2 08. Oxy- 
gen atoms are not shown for clarity. Grey dashed lines in- 
dicate chemical unit cell, (a) Polyhedral representation of 
structure. Light green and dark green polyhedra around Ba 
sites represent the two types of Ba coordination that alternate 
along the c-axis, 12 and 10-fold coordination respectively. Red 
tetrahedra illustrate the Mn'''^ coordination, (b) Simplified 
structure figure illustrating dominant exchange connectivity 
in Ba3Mn208 within and between the S = 1 dimers. Ji and 
J2 exchange correspond to distances of 4.569 and 5.711 A re- 
spectively. 



ing. Following the thermodynamic phase transition to 
lower temperatures, the experiments extrapolate to a 
zero-temperature ^iqHci = 9.04(0.15) T quantum criti- 
cal point. These measurements also indicate two phase 
transitions near the lower critical field implying the exis- 
tence of two magnetic long-range-ordered phases in close 
proximity to the quantum critical point. 

More recent measurements using single crystals con- 
firm the presence of two distinct ordered phases and 
provide a detailed map of the phase diagraroi^. EPR 
measurements indicate a zero-field splitting of the triplet 
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states, attributed in part to the efFects of single ion 
anisotrop}»i^, and implying more complex magnetic 
structures than previously suggested^^. 

The applicable Hamiltonian for this system is not yet 
fully characterized, and the magnitude and extent of 
dimer and interdimer exchange interactions are impor- 
tant for further understanding of Ba3Mn20g in both zero 
and applied magnetic fields. Here we describe tempera- 
ture dependent inelastic neutron scattering (INS) mea- 
surements on polycrystalline samples to determine the 
extent of dispersive magnetic excitations and exchange 
constants in Ba3Mn2 08- 

Ba3Mn208 is hexagonal (space group Rim) with 
room temperature lattice constants a = 5.711 and c = 
21.444 hM. The Mn5+ ions reside in a distorted tetra- 
hedral environment, as shown in Fig. [ija), resulting in 
an effective 5 = 1 moment. Dimer and interdimer 
magnetic interactions arc considered to be antiferromag- 
netic (J > 0) and Heisenberg with a small single-ion 
anisotropy, D. Including a Zeeman term for applied mag- 
netic fields H along the z-axis, the Hamiltonian is 

^ = E • S, +DY,{Stf g^^BHY^St, (1) 

i i 

where i and j designate coordinates of individual inter- 
acting spins, S. The nearest neighbor 5* = 1 sites along 
the c-axis, d ~ 3.985 A, have been proposed as strongly 
coupled antiferromagnetic S = I dimers. The dimers 
form an edge-sharing triangular lattice bilayer in the ab- 
plane. Neighboring planes along the c-axis are separated 
by alternating layers of oxygen coordinated Ba sites, ten 
and twelve site coordination respectively, as illustrated 
in Fig.[lja). Weaker interdimer exchange within the ac- 
plane, J2, and between the bilayers, Ji, have also been 
proposed based upon the crystal structure. These are 
illustrated in Fig. \V[h) along with the dimer exchange 
Jo. 



II. EXPERIMENTAL METHODS 

Powder samples of Ba3Mn2 08 were synthesized using 
stoichiometric amounts of BaC03 and Mn203 in a solid 
state reaction. The reactants were calcined under flowing 
oxygen at 900 °C for 30 h. The resulting green powder 
was then reground and sintered between 900 and 1000 °C 
under flowing oxygen for approximately 7 days. This 
growth procedure is similar to that described in Ref. [l^ . 
Single crystals of appropriate mass are unfortunately not 
yet available for studies of the dispersion using INS. 

SQUID magnetization measurements as a function of 
temperature for iioH = 1000 Oe and as a function of 
magnetic field for T = 40 K did not reveal any measur- 
able impurities due to Mn203 or Mn3 04 which are ferri- 
magnetic below Tc = 79 K and antiferromagnetic below 
Tc = 43 K respectivcl}^. The temperature dependent 
magnetic susceptibility compares well with previously 



published data with a rounded peak at T = 18 K and 
an activated low-temperature susceptibility characteris- 
tic of antiferromagnetic spin-gap systemsi^iii^. We also 
characterize the 150 < T < 350 K magnetic susceptibility 
via a Curie- Weiss law with Qqw = —43.2(2) K. 

Inelastic neutron scattering measurements were per- 
formed on 76 grams of Ba3Mn208 in an 18 mm diam- 
eter and 100 mm tall cylindrical aluminum sample can. 
The sample was produced from five separately prepared 
batches of Ba3Mn208. Each batch was checked for im- 
purity phases using SQUID magnetization and powder 
X-ray diffraction. Both measurements found no mea- 
surable impurity phases. Inelastic neutron scattering 
measurements were performed using the direct geome- 
try chopper spectrometer, DCS, at the NIST Center for 
Neutron Scattering. Spectra were measured for temper- 
atures between T = 1.4 and T = 160 K for one hour 
in each configuration. Temperature control was provided 
by a liquid He^ flow cryostat. Two incident wavelengths, 
A, were used. A = 2.9 A measurements probed energy 
transfers up to huj = 5.5 meV and wave vectors up to 
Q ~ 4.05 A~^ at the elastic position. A — 4.4 A mea- 
surements provide improved energy and wave vector res- 
olution up to Hlo = 3.05 meV and Q = 2.69 A~^. 
The calculated full width at half maximum (FWHM) 
energy resolution at the elastic position is Shio = 0.5 
and 0.15 fj.eV for the 2.9 and 4.4 A incident wavelengths 
respectively. Background measurements were made for 
each incident wavelength using an empty sample can at 
T = 1.4 K. A T < 200 K vanadium standard was mea- 
sured for calibration of detector sensitivity. A scattering 
angle dependent absorption correction for the cylindrical 
sample geometry was also applied to the inelastic scatter- 
ing intensity. Unless otherwise noted, these backgrounds 
and normalization are applied to our presented results. 

Neglecting Bose occupation and Debye- Waller factors, 
the magnetic neutron-scattering cross section is propor- 
tional to the scattering function, iS(Q,ti;), 

^4|F(Q)P5(Q,.), (2) 

where k' and k are the magnitude of the final and initial 
neutron wave vectors and F{Q) is the magnetic ion form 
factor. We plot our measured scattering intensity in units 
of iS(Q,ct;), but we do not normalize the data by the 
magnetic form factor. 

III. RESULTS AND DISCUSSION 

Figure [2] shows the scattering intensity as a function 
of wave vector transfer in the vicinity of the elastic posi- 
tion, —0.1 < TiLo < 0.1 meV, at several temperatures. 
Nuclear Bragg peak positions and intensity compare 
well with the previously determined room-temperature 
structurei^. We note that for T < 80 K there are two ad- 
ditional weak Bragg peaks at Q = 1.096 and 1.545 A~^. 
These arc likely associated with either Mn2 03 or Mn304 
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FIG. 2: (color online) Temperature dependent diffraction for 
Ba3Mn208 measured in the A = 2.9 A configuration. Data are 
integrated between -0.1 and 0.1 meV energy transfer. Data 
have not been background subtracted or corrected for absorp- 
tion. Higher temperature data are offset vertically from the 
T = 1.4 K data for presentation. Several characteristic Bragg 
peaks are indexed in the figure. 



impurities as discussed earlier. We do not observe any 
contribution from spin-waves due to these impurities in 
the inelastic portion of the spectra. 

Although wave vector resolution was not optimized for 
diffraction w 0.03 for A = 4.4 A), we fit the (003), 

Q=1A 

(006), (009) and (110) Bragg peak positions for each inci- 
dent wavelength to determine lattice constants as a func- 
tion of temperature, c.f. Fig. [3^a) and (b). These val- 
ues are consistent with the previously determined room 
temperature structure. There are no apparent structural 
phase transitions from T = 160 to 1.4 K. Both the a and c 
lattice vectors contract at lower temperatures, but there 
is only a 0.3% change in the a-axis lattice constant and an 
even smaller change in the c-axis lattice constant, 0.06%, 
from T = 295 K to T = 1.4 K. Linear fits to these data 
along with the T = 295 K values are shown in Fig. [31^ a) 
and (b) along with calculated values of c/a and the re- 
spective calculated curve in panel (c). The fitted lines 
provide a good description of the data with coefficients 
of linear expansion: Uc ~ c(t^o) IT ^ 2.3(l)i? — 6 and 
- - 1 ''" ^ 1.25(1)^;- 5 K-i. Including quadratic 



"1 ~ a(T=0) dT 

terms does not substantially improve the fits. The weak 
lattice parameter temperature dependence indicates that 
any changes in the magnetic excitation spectra with tem- 
perature are not likely associated with changes in ex- 
change due to changes in distance between individual 
spins. 

The T = 1.4 K excitation spectra for the A = 4.4 
and 2.9 A configurations arc shown in Fig. HJa) and (b). 
There is a single dispersive excitation with a gap of A « 
1 meV and a bandwidth of approximately 1.5 meV. The 
decreasing scattering intensity with increasing wave vec- 
tor transfer immediately suggests the spectrum is mag- 
netic. The increase in scattering intensity near Tiuj = 0.25 



21.47 
--21.45 
o 21.43 

21.41 
5.71 

< 5.70 

CD 

5.69 
3.77 
" 3.76 
3.75 



Ref. 18 



H 1 1 1 1 1 1 W 




^ — I — 1 — I — 1 — I — 1 — W NHH 




40 



80 
T(K) 



120 



16 



FIG. 3: (color online) Temperature dependent lattice param- 
eters of Ba3Mn20g. (a) c-axis lattice constant versus tem- 
perature, (b) o-axis lattice constant versus temperature, (c) 
c/a as a function of temperature. Solid points from Ref. [l3 
for T = 295 K are shown for comparison. Lattice constants 
were obtained from the weighted mean peak position of the 
(003), (006), (009) and (110) Bragg peaks for A = 2.9 and 
4.4 A. Error bars are the error in the mean. Solid lines in (a) 
and (b) are linear fits described in the text. Solid line in (c) 
is calculated based upon the linear fits for the a and c-axis 
lattice constant. 



meV in Fig. [H^a) for small wave vectors is not intrinsic 
to the sample and has been seen in background measure- 
ments. Figure and (f) show the low-temperature 
wave vector integrated scattering intensity. There is a 
peak in the magnetic density of states in the vicinity 
of the spin-gap energy, but there is no peak at the top 
of the dispersive band of excitations. The spectrum of 
one-dimensional Heisenberg gapped quantum antiferro- 
magnets has characteristic Van-Hove singularities at the 
top and bottom of the band. The absence of any peak 
at higher energy transfers implies the magnetic excita- 
tions in Ba3Mn208 are at least two-dimensional in their 
connectivity. 

For comparison to our results, we calculate the INS 
cross section of = 1 antiferromagnetic dimers with 
weak interdimer interactions. At finite temperature the 
scattering function is written as2£ 



Z 
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(3) 



where f3 = jt^, Z is the partition function and the sum 
is over the initial and final states ip and tp' with energies 
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FIG. 4: (Color online) Measured and calculated INS intensity for Ba3Mn208 at T = 1.4 K. Panels (a) and (b) are data 
measured for incident wavelengths of A = 4.4 and 2.9 A respectively. Panels (c) and (d) are the calculated scattering intensities 
corresponding to the coupled S = 1 dimer model described in the text. Panels (e) and (f) are the energy dependent integrated 
scattering intensity of the data and model calculations shown in panels (a)-(d) for wave vectors 0.65 < Q < 2 and 
0.75 < Q < 3.25 for A = 4.4 and 2.9 A respectively. Model lineshapes in (e) and (f) are fit to the data including an elastic 
Gaussian lineshape and a constant background as discussed in the text. 



and E^i. For a dimcr, Eq. [3] becomes 



5e'^^[l-cos(Q.d)] ^^^ 



5{E^ - E^, + fiuj), 



(4) 



where xi and X2 are the respective crystallographic co- 
ordinates of the spins in the dimer and Sn are the spin 
operators. 

An isolated antiferromagnetic S" = 1 dimer with in- 
tradimer exchange Jq will have a non-magnetic (total 
spin St = 0) singlet ground state at an energy of 
E = —2 Jo with St = ^ triplet and St = '2 quintet excited 
states at energies E = —Jq and E = Jq. This results in 
triplet and quintet spin-gaps of hio = Jq and huj = 3,/o- 
However, magnetic INS will only probe the singlet-triplet 
and triplet-quintet cross-sections, i.e. \ASt\ = 1. In- 
cluding structure factors from Eq. IH the scattering func- 
tion is a sum of two terms, 



where d is the bond vector between the spins of the 
dimer. Although the calculated matrix elements of the 
triplet-quintet transitions are larger than the singlet- 
triplet transitions, the triplet-quintet transitions are 
thermally activated with less spectral weight and will 
only be populated at higher temperatures. 

Equation [5] does not account for interdimer correla- 
tions, i.e. dispersive excitations. The random phase 
approximation (RPA) has been successful in describing 
the dispersion of weakly coupled dimers in spin-gap sys- 
tems. This has been illustrated for varying numbers of 
interactions and spin-quanta in several experimental sys- 
tems including KCuClg^ii^, TlCuCl^^, CsgCraBrc^, 
CsaCrzl!^, PHCCi and BaCuSiaOg^. The RPA dis- 
persion for Heisenberg exchange coupled dimers is 



hiviQ) = + A/2A J(Q)i?(T) 



(6) 



5(Q,c^) 



4e2'''/^[l-cos(Q-d)] 



S{huj -J)+ (5) 



where M is the transition matrix element (M — 
1]), .^(Q) is the Fourier sum over interactions beyond 



\S[S~ 
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dimcr exchange, A = Jg and R{T) is the thermal popu- 
lation difference between the ground and excited states. 
For S = 1 antiferromagnetic dimers, 



RiT) 



1 - exp(-A/?) 



1 + 3 exp(- A/3) + 5 exp(- A/3) ' 



(7) 



considering only singlet-triplet excitations. Ba3Mn2 08 
has a single dimer per unit cell and the interdimer in- 
teractions propagate the triplet excitation leading to the 
Fourier sum 



LU2 ~ cos(27rfc) -t- cos(27r[/i + fc]) + cos(27r/i) 
LOi = cos(— [— /i -I- fc -I- Z]) -I- cos(— [— /i — 2fc - 



27r 

-|-cos( — [2h + k - 
3 

J{Q) = 2J2CJ2 + Jiwi. 



I]) 



(8) 



Recent EPR measurements have revealed a zero-field 
splitting of £) = —0.032 meV, although both modes will 
have an identical dispersioni^. Examination of thermo- 
dynamic measurements have included an exchange con- 
stant, J3, which couples neighboring dimers in a bilayer 
from spin-1 of a dimer to spin- 2 of a second dimer—. 
This exchange term represents a change in phase of the 
triplet excitation between dimers in the ab plane, and 
would change the prefactor of the UJ2 term in Eq. [8] to 
be 2(J2 — J3). We have chosen to omit the J3 exchange 
from the current analysis. Because of its large spin-spin 
distance (6.964 A) and out of plane coupling, it is pre- 
sumably much weaker than J2 and Ji . In addition, recent 
calculations examining the relative strength of exchange 
constants in Ba3Mn208 have shown that the J3 exchange 
constant is effectively zero^^. 

The scattering function must also be modified to 
account for the dispersive excitations. The single 
mode approximation (SMA) has been successfully ap- 
plied to several dispersive disordered gapped antifer- 
romagnets with interdimer exchange included in the 
Hamiltonian^>^i21i^i^. This results in an additional 
multiplicative ^^^qj term in 5(Q,w) such that consid- 
ering only singlet-triplet scattering for Ba3Mn208 the 
scattering function becomes 



5(Q,c^) = 



4e2-^'«[l -cos(Q -d)] 



-5{fiuj-fiuj{q,)), 



(9) 

where 7iw(Q) is given by Eqs. [6]l8l The SMA is appro- 
priate for the case of Ba3Mn208 for T ^ Jq since only 
singlet-triplet excitations will be effectively probed via 
INS in this temperature regime. 

For comparison of our high resolution (A = 4.4 A) 
polycrystalline measurements we numerically spherically 
average Eq. [21 



5(Q,a;)= / ^^(Q,^). 



(10) 



This process was also recently employed successfully in 
examination of polycrstyalline measurements of a gapped 
antiferromagnet using the same instrumentation^!. The 
interpretation of our measurements is more straightfor- 
ward given the absence of hydrogen or contamination 
from phonons in our spectrum. We calculate 5(Q,a;) 
over spherical shells in |Q| space at fixed values of en- 
ergy transfer with d fixed as the proposed dimcr bond 
vector. This was done for a series of Jq, Ji and J2 values. 
This spectrum was combined with an identical spectrum 
shifted in energy transfer by the value D = —0.032 meV. 
The spectrum was then multiplied by \F{Q)\'^ g^^^ 
convolved with a Gaussian representation of the mean 
instrumental energy and wave vector resolution over the 
energy and wave vector range of the magnetic excitation. 
A constant background and multiplicative prefactor were 
used as fitting parameters of the calculated spectrum 
in comparison to the measured data. This procedure 
yields best fit exchange constants Jq — 1.61(3), Ji — 
~0M2toml and J2 = O.lUtoml meV^. The corre- 
sponding best fit INS scattering intensity for both instru- 
ment configurations is plotted in Fig. \M,c) and (d) , and 
agrees very well with the dispersion and intensity mod- 
ulation observed in the measurement. The determined 
value of Jo and the corresponding spin-gap based upon 
the dispersion, A = 1.05 meV, are both in the vicinity of 
values from thermodynamic measurements. 

Figures IDJe) and (f ) show the fitted wave vector inte- 
grated lineshapes compared to the measured data. Fit 
parameters for the calculated lineshapes include an elas- 
tic Gaussian, an overall multiplicative prefactor and a 
constant background. The calculations based upon both 
incident wave lengths agree very well with the measure- 
ment. 

We plot the low-temperature singlet-triplet dispersion 
relation in Fig. [5l Overall minima in the dispersion occur 
at the (| — (5, | — ^, and (| + 15, | + wave vectors 
for I = 3n and (^ + 5, ^ + S, I) and (| - (5, | - (5, /) 
wave vectors for I — | + 3n where n is an integer and 
6 ~ 0.025. The nonzero value of Ji results in a finite 
dispersion along the {(^ ^ I) direction with a periodicity 
of three reciprocal lattice units (rlu). This is shown in 
Fig. \5lh) for the curves plotted using the top axis. 

The powder average of the first frequency moment of 
the T = energy integrated scattering function. 



h{uj)Q = 



dnduj 



An 



nuj{ci)s{ci,nLu) 



« iF(Q)ni 



sm{Qd) 
Qd ^ 



(11) 



provides direct information regarding the length of the 
dimer bond, d. Figure [B] shows the first moment of the 
measured T = 1.4 K scattering intensity for both incident 
wavelengths. These data are fit to Eq. [11] including an 
overall constant and multiplicative prefactor. A simulta- 
neous fit of both data yields a good representation of the 
measured results with d = 4.073(7) A. This value agrees 
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FIG. 5: (Color online) T = 1.4 K singlet-triplet single crys- 
tal dispersion relation determined for Ba3Mn2 08 based upon 
powder average analysis, (a) Contour figure of dispersion in 
the hexagonal (HKO) plane, (b) Dispersion as a function of 
reciprocal lattice units (rlu) along other primary directions of 
the hexagonal structure. Dispersion along the (0 L) and 
(1 L) directions is plotted using the top axis. 




FIG. 6: (Color online) First moment of scattering intensity 
integrated between Tiu) = 0.75 and 3 meV for the T = 1.4 K 
A = 4.4 and 2.9 A measurements shown in Fig. |4ja) and (b). 
Solid lines correspond to fits described in the text with the 
fitted |d| value of d = 4.073(7) A. 



with the description of the dimer bond being the short 
vertical bond between Mn^+ moments in the Ba3Mn2 08 
crystal structure d = 3.985 A, c.f. Fig. [TJ 

The magnetic spectra of Ba3Mn208 as a function of 
temperature for 1.4 < T < 160 K was also measured. In 
order to consider a range of energy transfers beyond the 
singlet-triplet band, we examine the A = 2.9 A configu- 
ration data. Constant wave vector scans for 0.9 < Q < 
1.0 A-i are shown in Fig. [7| for several temperatures for 
both neutron energy gain and loss. This wave vector 
was chosen because it is in the vicinity of both the peak 
in the density of states and the overall minimum in the 
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FIG. 7: (Color online) Energy dependent scattering inten- 
sity for 0.9 < Q < 1.0 A^^ for a series of temperatures as 
measured in the A = 2.9 A configuration. Data for differ- 
ent temperatures are vertically offset. Solid lines are fits to 
the Lorentzian functions described convolved with the energy 
dependent resolution of the instrument. Inset illustrates scat- 
tering intensity of the T = 12 K measurement plotted on a 
logarithmic intensity scale. Dotted red line corresponds to 
the spectral contribution of the lower energy mode and the 
dashed blue line corresponds to the higher energy mode. 



measured dispersion. From the T = 1.4 K spectrum in 
Fig. m^b) , this wave vector also provides minimal inter- 
ference from higher energy portions of the singlet-triplet 
spectrum. The T ~ 1.4 K spectrum consists of a single 
peak at htu « 1.25 meV with no significant scattering in- 
tensity on the neutron energy gain side of the spectrum. 
As temperature is increased, the single peak broadens, 
moves to slightly larger energy transfers and additional 
scattering intensity develops in the vicinity of 3 meV. 
There is also additional scattering intensity which devel- 
ops on the neutron energy gain side of the spectrum at 
elevated temperatures. 

In order to distinguish singlet-triplet and thermally ac- 
tivated triplet-quintet excitations, we model the constant 
wave vector scans using two excitations. We fit to two 
Lorentzian functions of the form^ 



(«M + i>( 



AT, 



AFi 



SFa 



(12) 



where Fi, wi, F2 and lu2 a-i'c the half width at half 
maximum and energy of two respective excitations with 
the Lorentzian width reflecting temperature dependent 
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FIG. 8; (Color online) Temperature dependent integrated in- 
tensity and characteristic energy determined from fits to data 
shown in Fig. [7] (a) Total and neutron energy gain and 
loss portions of the integrated scattering intensity for both 
modes. The solid black line is a fit of the temperature depen- 
dent total scattering intensity of the isolated antiferromag- 
netic S = 1 dimer including both singlet-triplet and thermally 
activated triplet-quintet excitations. The dashed line only in- 
cludes singlet-triplet excitations. Solid blue and green curves 
are calculated from the singlet-triplet and triplet-quintet fit- 
ted scattering intensity, (b) Mode energy as a function of 
temperature. The T — 1.4 K spin-gap value determined from 
fitted exchange parameters is shown as a single red diamond. 
Solid line is a calculation of the RPA model temperature de- 
pendent spin-gap using the determined values of exchange 
constants. 

broadening of the spectrum. The Bose factor {n{Lu) + l) = 
[1 — exp(— enforces detailed balance of the scat- 
tering intensity. Equation [T^] and a variable width and 
amplitude Gaussian peak at the elastic position were con- 
volved with the energy transfer dependent instrumen- 
tal energy resolution with fitting parameters A, B, Ti^2 
and A time independent background of the form 

BG oc (^jT^ — hu!)~'^, as typically used for time-of-flight 
direct geometry INS measurements and a constant back- 
ground were determined from the T < 40 K data and held 
fixed. The two excitation fits are shown as solid lines in 
Fig. [T] The lineshapes agree well with the energy depen- 
dent distribution of scattering intensity for both neutron 
energy loss and gain. To further illustrate the presence 
of the triplet-quintet excitation, we plot the T = 12 K 
data and fitted lineshape including the individual contri- 
butions from each mode in the inset of Fig. [T] 

Figure[S]^a) shows the total integrated scattering inten- 
sity for both modes as well as the neutron energy gain 
and loss portions as a function of temperature. We fit 



the total scattering intensity to the temperature depen- 
dent isolated dimer scattering function in Eq. O with an 
overall multiplicative prefactor and the exchange Jo as 
fitting parameters. Considering only singlet-triplet exci- 
tations yields the dashed red line in Fig.tSja) and a value 
of Jo = 3.9(2) mcV. However, including both singlet- 
triplet and triplet-quintet modes yields a much better 
fit (solid black line) and a more appropriate value of 
Jo = 1.9(3) mcV. The improved fit implies that there is a 
non-negligible contribution of thermally activated triplet- 
quintet excitations at higher energy transfers. The calcu- 
lated temperature dependent scattering intensity for the 
neutron energy gain and loss portions of the spectrum 
based upon the singlet-triplet, triplet-quintet model also 
agrees well the data shown in Fig.[S]Ja). 

The fitted values uji and W2 are plotted in Fig. [5]^b). 
The LOi value is not the spin-gap energy, rather it is 
the characteristic energy of the lower energy mode found 
for the constant wave vector scan. For comparison, we 
also plot the determined spin-gap energy based upon the 
T = 1.4 K powder average analysis. The calculated tem- 
perature dependent spin-gap of the singlet-triplet exci- 
tation based upon Eqs. [5][5] and the fitted exchange con- 
stants is plotted as a solid line in Fig. [HUb). The RPA 
coupled dimer description agrees well with the tempera- 
ture dependence of the energy scale of the singlet-triplet 
excitation. We also note that above base temperature the 
RPA function agrees with the temperature dependence 
of the activated triplet-quintet excitation, lu2, although 
shifted to higher energy transfers, implying that these 
excitations may share a similar temperature dependent 
dispersion renormalization. 



IV. CONCLUSIONS 

Through INS measurements we have shown that 
there exists a well-defined singlet-triplet spectrum in 
Ba3Mn2 08. Although the measured bandwidth of the 
magnetic spectrum is larger than the spin-gap, the ex- 
change constants indicate that Ba3Mn208 can be consid- 
ered a triangular lattice of weakly coupled 5" = 1 dimers. 
Comparison to an appropriate RPA coupled dimer de- 
scription of the scattering function is able to determine 
the exchange constants. In addition, the wave vector 
dependent scattering intensity agrees with the dominant 
dimer bond being the predicted short vertical bond illus- 
trated in Fig. [Ijb). 

The examination of temperature dependent scattering 
intensity indicates that both singlet-triplet and triplet- 
quintet excitations are observed in the INS spectrum 
(there is no INS cross-section for S* = 1 antiferromag- 
netic singlet-quintet excitations). Based upon the rel- 
ative energy scales of the two observed excitations, we 
can estimate the mean energy of singlet-quintet excita- 
tions as w 3.4(4) meV. Single-crystal INS measurements 
may be able to determine the dispersion associated with 
thermally activated triplet-quintet excitations or perhaps 
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observe multi-particle excitations, quintet-triplet decay 
or interference of single- and multi-particle excitations^^. 
The existence of singlet-quintet and triplet-quintet ex- 
citations may also be able to explain the heat capacity 
above T « 3 K which can not be accounted for by single 
excitation models^^. 

We also point out that the currently determined 
exchange constants may be able to further describe 
the magnetic field dependent phase diagram at low- 
temperatures or place limits on the nature of the pro- 
posed long-range-ordered phases. The difference between 
the observed phase diagrams for H\\c and H J- c must be 
related to the single ion anisotropy term in the Hamilto- 
nian, Eq. [TJ As we will show elsewhere^^, this term in- 
duces an effective exchange anisotropy in the low-energy 
Hamiltonian that results from projecting the original 
Hamiltonian into the subspace generated by the singlet 
and Sz = 1 triplet of each dimer. These two states can 



be described with a pseudospin 1/2 variable. The com- 
bined effect of geometric frustration and anisotropy leads 
to the appearance of a new phase for _L c that will be 
discussed in Ref. [13. 
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